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SHAMO:  A  COMPUTER  PROGRAM  FOR  ESTIMATING  THE  POSITION  AND  SHAPE 
OF  THE  SURFACE  SHADOW  ZONE  IN  SONAR  OPERATIONS 

by 

Bernard  de  Raigniac  an'’.  John  Padley 


ABSTRACT 


SHADZO  is  a  composite  program  which  calculates  the  position  and 
shape  of  the  first  surface  shadow  zone  as  a  function  of  the  source 
depth  for  a  given  sound  speed  profile.  The  program  is  small 
enough  to  be  run  on  a  shipboard  mini-computer  in  a  few  minutes. 


I 


INTRODUCTION 


It  has  been  found  from  SACLANTCEN’s  reverberation  studies  that 
information  on  the  position  and  extent  of  the  surface  shadow  zone* 
is  needed  during  at-sea  experiments.  Thus,  a  program  has  been 
developed  which,  given  the  sound  speed  profile,  calculates  the 
pertinent  characteristics  of  the  shadow  zone  and  can  be  run  on  a 
shipboard  mini-computer  in  a  few  minutes.  In  orde  '  to  make  the 
maximum  possible  information  available  to  the  scientists  aboard 
ship,  the  program  also  includes  the  estimation  of  shadow  zone 
shape. 

SHADZO  (surface  SHADow  ZOne  program)  is  a  composite  program  which 
calculates  the  position  and  shape  of  the  surface  shadow  zone, 
as  a  function  of  the  source  depth,  for  a  given  sound  speed  profile. 
These  characteristics  are  obtained  by  direct  investigation  of  the 
limiting  rays,  rather  than  by  representation  of  the  sound  field 
by  ray-tracing  with  a  high  density  of  rays. 

This  memorandum  describes  the  theoretical  basis  for,  and  the 
implementation  of,  the  two  basic  parts  of  SHAZDO:  (1)  a  part  that 
calculates  the  distance  (i.e.,  the  inner  range  limit),  the  extent 
and  the  maximum  thickness  as  functions  of  source  depth,  and  (2)  a 
part  that  estimates  the  shape  of  the  shadow  zone  for  selected  source 
depths.  A  flow  chart,  a  listing,  and  explanatory  diagrams  of  SHAZDO 
are  given  in  Appendix  A. 


The  surface  shadow  zone  treated  here  is  the  first  (i.e,,  shDrtest-range 
shadow  zone;  the  recurring  surface  shadow  zones  a c  longer  ranges  are  not 
covered. 


1. 


THEORY:  POSITION  AND  SHAPE  OF  THE  SURFACE  SHADOW  ZONE 


Background 


When  the  sound  source  is  located  above  the  critical  depth",  zCritJ 
there  is  a  shadow  zone  along  the  surface  at  certain  ranges.  This 
surface  shadow  zone  has  been  discussed  by  Mellberg,  who  calls  it 
"shadow  layer",  in  Ref.  1  .  Mellberg  provides  formulae  for  the 
maximum  thickness  and  horizontal  extent  of  the  shadow  zone,  based 
cn  a  2-layer,  constant  gradient  model  of  the  medium.  The  present 
memorandum  provides  a  further  description  of  the  characteristics 
of  the  shadow  zone,  based  on  a  multi-layer,  constant-gradient 
model . 


Figures  1  and  2  illustrate  the  formation  of  the  surface  shadow 
zone  in  a  medium  approximated  by  a  two-layer,  constant-gradient 
model.  Figure  la  shows  the  rays  which  limit  the  shadow  zone, 
both  in  range  and  depth,  when  the  source  is  between  the  surf 
and  the  minimum- speed-depth,  zm<  At  the  surface,  the  shadow  zone 
is  bounded  by  rays  which  have  zero  grazing  angle.  The  maximum 
depth,  or  maximum  thickness  of  the  shadow  zone,  is  the  depth  at 
which  the  ray  that  was  horizontal  at  the  source  becomes  horizontal 
again.  When  the  source  is  located  between  the  mi nimum-speed  depth 
z  and  the  critical  depth,  z  ...  as  in  Fig.  lb,  a  shadow  zone 
of  smaller  range  extent  and  smaller  maximum  thickness  will  occur. 


Figure  2  illustrates  the  distance,  extent  and  maximum  thickness  of 
the  shadow  zone,  again  for  a  2-layer  model.  The  multi-layer  model 
will  be  introduced  next,  and  then  these  three  characteristics  will 
be  discussed  in  turn. 


T^e  critical  depth  is  the  aepth  at  which  the  soond  speed  is  the  same 
as  at  r<~e  sea  s^-r^ace 


In  SHADZO,  ?.  multi-layer,  constant-gradient  model  of  the  medium 
is  used.  The  rays  propagate  along  circle  segments  in  each  layer, 
as  illustrated  in  Fig.  3*  The  radius  of  the  path  in  the  ith 
layer  is  given  by: 

R  ■=  Jl 

i  S±  CEq-  !] 

where 

k  =  Snell’s  constant  of  the  ray 
=  c/cosft  (at  any  depth) 

g.  =  gradient  in  the  ith  layer 


The  horizontal  distance  travelled  by  the  ray  in  the  ith  layer  is 
then: 

Ad±  -  Ri(sinai_1  -sinc^) 


[Eq.  2] 


1 . 2  Distance 

The  distance  to  the  shadow  zone-.  >s  illustrated  in  Fig.  2,  is  defined 
as  the  inner  range  limit  of  the  shadow  zone  at  the  surface. 

Appendix  A  of  Ref.  2  provides  formulae  for  the  distance,  for  both 
2-Iayer  and  r..ulti-layer  constant-gradient  models. 


From  Eq.  2,  we  see  that  the  total  horizontal  distance,  d,  travelled 


by  a  ray  in  traversing  n 
the  surface  is  given  by 


n  n 

=  S  fid.  =  E 
i=l  1  i~l 


Azi 


layers  between  the  source  depth  and 


-  C. 

1-1 


d 


c . 
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To  avoid  computing  errors  when  c^sac.  this  expression  may  be 
transformed  into 


[Eq.  3] 


d  may  be  expressed  as  a  function  of  grazing  angle  at  the  surface,  y, 
by  writing  k  =  c0/cos  y  -  Thus,  we  see  that  the  distance  to  the 
shadow  zone,  D,  is  obtained  from  Eq.  3  simply  by  setting  Y  =  0» 
then,  we  have  k  =  c0  and 


D  = 


n 

2 


i=l 


[Eq.  4] 


1 . 3  Extent 


The  horizontal  extent  of  the  shadow  zone  at  the  surface  is  determined 
by  two  rays  having  zero  grazing  angle,  as  illustrated  in  Fig.  2. 

It  can  be  seen  from  Fig.  2  that  the  horizontal  distance  travelled 
by  the  long-range  limiting  ray  in  going  between  source  depth  and 
critical  depth  is  a  measure  of  the  extent  of  the  shadow  zone. 

Thus,  the  extent  is  given  by 


E 


k 

2  2 
i=n+l 


[Eq.  5] 


where  it  is  assumed  that  there  are  n  layers  between  the  source 
depth  and  the  surface  and  k  layers  between  the  critical  depth 
and  the  surface. 


. 


1.4  Maximum  Thickness 

As  discussed  ip  Ref.  1  and  illustrated  in  Figs.  1  and  2,  the 
maximum  depth  of  the  shadow  zone  is  the  depth  where  a  rty 
that  was  horizontal  at  the  source  becomes  horizontal  again. 

At  this  depth,  z.  ,  the  sound  speed  is  equal  to  c  ,  the  sound 
speed  an  source  depth. 

1. 5  Shape 

The  shape  of  the  shadow  zone  can  be  defined  as  the  envelope  of 
the  downward- refracted  rays,  as  illustrated  by  the  dashed  curves 
in  Fig.  1.  This  envelope  is  relatively  difficult  to  calculate 
exactly,  so  an  estimate  of  the  shape  is  obtained  as  follows. 
First,  at  close.y- spaced  depths,  we  determine  the  lcrus  of 
points  at  which  r ’.ys  vertex  (i.e.  become  horizontal).  For 
the  case  shown  in  Fig.  lb,  we  see  that  this  locus  would  provide 
a  reasonable  estimate  of  the  envelope *.  Figure  4  illustrates 
the  case  of  a  limiting  ray,  such  as  might  be  encountered  in 
the  presence  of  a  strong  negative  gradient;  in  this  case,  the 
locus  of  vertices  is  not  a  good  estimate.  To  improve  the 
estimate,  at  each  depth,  the  intersections  of  rays  vertexing 
at  shallower  depths  are  determined.  Finally,  at  each  depth, 
that  point  (vertex  or  intersection)  is  used  which  resuits  in  the 
smallest  extent  of  the  shadow  zone.  The  estimation  of  sho^  w 
zone  shape  from  these  vertices  and  intersections  is  illustrated 
in  Fig.  A. 


It  would  be  clearly  not  provide  a  reasonable  estimate  for  the  case  in  Fig.  la; 

nowever,  this  case  is  not  addressed  here,  as  the  snaps  estimation  part  of 

SHADZO  nas  been  implemented  only  for  cases  in  which  z  >z  , 

s  m 
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FIG.  2  DISTANCE,  EXTEN  f  AND  MAXIMUM  THICKNESS  OF  THE  SHADOW  ZONE 


2. 


DESCRIPTION  OF  THE  SHADZO  PROGRAM 


2.1  Hardware  Configuration 

The  program  was  written  to  run  under  the  Hewlett-Packard  Real-Time 
System  [Ref.  3]>  and  thus  requires  a  minimum  hardware  configuration 
of : 


H-P  2116B  Computer  with  16K  Memory 
H-P  12578A  Direct  Memory  Access 
H-P  12579A  Extended  Arithmetic  Unit 
H-P  12591A  Memory  Protect 
Fixed  Head  Disc  or  Drum  Storage  Unit 
Time  Base  Generator 

Operator  Console  (ASR-33  or  ASR-35  Teleprinter) 

In  addition,  the  SHADZO  program  requires  a  Tektronix  T4002/4802 
Graphic  Computer  Terminal  with  Tektronix  Joystick  015-0175000  and 
Hard  Copy  Unit  46OI .  The  Tektronix  Terminal  can  be  used  in  place 
of  the  operator  console,  if  an  additional  paper-*tape  input  device 
is  available.  In  the  absence  of  a  Tektronix  Terminal,  the  program 
can  still  produce  a  printed  output  on  an>  output  writer. 

2 . 2  Sottware  Configuration 

The  program  requires  a  Real-Time  System  generated  to  provide  a 
minimum  background  disc-resident  area  of  20K  (octal)  locations. 

The  Tektronix  Terminal  should  be  allocated  logical  unit  number  16 
and  should  be  used  with  the  Real-Time  teleprinter  drivei  DRV00 . 


The  data  are  input  to  the  program  via  a  punched  tape  containing 
the  temperature  or  sound  speed  profile.  Table  1  gives  an  example 
of  such  a  profile,  as  well  as  the  required  format. 

TABLE  1 

TYPICAL  DATA  TAPE 


FL  11- - 

29,7,71,1,40 

_16  - - 

0  ,24.4 
6  ,23.5 
8  ,22.8 
12, 22. 2 

16.19.5 
18,18  **“ 

20.16.5 

24.15.5 
28,14.8 
40,14.3 

64.13.5 

180.13- 5 

250.13- 45 

400.13- 4 

450.13- 35 
2500,13.35 


Profile  name— up  to  6  alphanumeric  charact* 

•Day,  Month,  Year,  Hour,  Minute—  5  integers. 

No.  of  points  on  the  profile  -  integer, 
negative  for  temperatures,  positive  for 
.velocities. 

Temperature  or  velocity  profile. 

Pairs  of  points  giving  depth  in  metres  and 
temperature  in  °C  or  velocity  in  m/s  . 


In  the  conversion  from  temperature  profile  to  sound  speed  profile 
[according  to  the  Leroy  formula  (Ref.  4)],  the  salinity  and 
longitude  are  considered  constants,  being  pre-set  to  38.6  parts/ 
thousand  and  40°,  respectively.  Source  depth  values  are  input  via 
the  Operator  Console. 


2.4  Program  Details 

The  SHADZO"  program  is  trritten  in  Fortran  II,  with  the  exception 
of  the  routines  for  the  Tektronix  Terunal,  which  are  standard 
SACIAKTCES  library  routines  written  in  Assembler  Code. 

The  program  consists  of  four  basxc  sections,  as  follows: 

s 

Reading  of  the  data  tape  t ape,  ar  d  conversion 
of  temperatures  to  sound  speeds  if  necessary. 

Calculation' and  printing  of  distance,  extent 
and  maximum  thickness. 

Calculation  and  graphical  display  of  distance, 
extent  and  maximum  thickness. 

Calculation  and  graphical  display  of  shadow 
zone  shape . 

The  first  section  is  performed  once  for  each  set  of  data;  the 
other  sections  are  performed  as  many  times  as  required. 

Sections  two  and  three  can  be  run  with  two  alternative  models 
of  the  sound  speed  profile:  (1)  a  multi-layer,  constant-gradient 
model,  using  all  of  the  sound  speeds  calculated  from  the  Leroy 
formula;  and  (2)  an  approximation  using  two  constant-gradient 
layers . 


Section  One 


Section  Two: 


Section  Three: 


Section  Four: 


Because  of  the  restriction  m  H-P  Fortran  II  which  allows  only  five  characters 
per  identifier,  the  program  name  actually  used  is  "SHAZO"  .  There  are  some 
discrepancies  between  the  terminology  used  in  the  text  of  this  memorandum  and 
the  terminology  usea  in  the  program  and  its  outputs  shown  below.  The 
correspondence  is  as  follows  : 

Term  in  Text 

Souna  speed 
Maximum  thicxress 
Distance 
Extent 


Term  3n  Program 

Velocity 

Thickness 

Range 

Extension 
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A  flow  chart  and  listing  of  the  program,  together  with  diagrams 
indicating  the  significance  of  variable  names  for  points  on  the 
s*Hind  speed  profile,  are  given  in  the  appendix. 


The  outputs  from  the  program  are  prints  and  plots,  as  described 
below. 

2.5.1  Prints 

a.  Sound  speed  profile:  A  sample  output  is  given  in 
Table  2v 

TABLE  2 

EXAMPLE  OF  SOUND  SPEED  PROFILE  OUTPUT 


SHADOW'  ZONE  PREDICTIONS 


XBT  FL  11 

DATE  29-7-71 

TIME  1-40 

DEPTHS 

VELOCITIES 

METRES 

METRES/SEC 

.00 

1557-3 

6.00 

1535.2 

8.00 

1533-5 

12.00 

1532.1 

16.00 

1525.0 

18.00 

1520.8 

20.00 

1516.4 

24.00 

1513.5 

28.00 

1511.4 

40.00 

1510.0 

64.00 

1507.8 

180.00 

1509.7 

250.00 

1510.7 

400.00 

1513.0 

450.00 

1513.7 

2500.00 

1548.0 

b.  Distance,  extent-  and  thickness:  An  example  is 

given  in  Table  3  for  the  multi-layer  model.  The  output  froa  the 
two-layer  aodel  has  a  sinilar  fonsat. 


TABLE  3 

EXAMPLE  OF  SHADOW  20KE  CHARACTERISE  ICS  OUTPUT 


MULTILAYER  MODEL 


SOURCE  DEPTH 

EXTEKSIOK 

RASGE 

THICKS ESS 

S.  VELOCITY 

METRES 

METRES 

METRES 

METRES 

M/SEC 

100.0 

35826.4 

794-3 

57-5 

1508.40 

20«?.0 

34794-9 

1310.0 

39.8 

1510-00 

300.0 

33733-8 

1840.6 

27.8 

1511.48 

400.0 

32640.1 

2387-4 

24-9 

1513-02 

500.0 

31511-9 

2951-6 

22.6 

1514-52 

600.0 

30341-5 

3536.8 

20.3 

1516.20 

700.0 

29122.7 

4146.2 

19-3 

1517.87 

-00.0 

27849-2 

4782.9 

18.6 

1519-55 

900.0 

26513-1 

5451-0 

17-8 

1521.22 

1000.0 

25104-4 

6155.3 

17-0 

1522.90 

1100.0 

23610.1 

6902.5 

16.2 

1524.57 

1200.0 

22012.7 

7701.2 

15.3 

1526.25 

1300.0 

20288.1 

8563.5 

14-4 

1527.92 

1400-0 

18400-3 

9507.4 

13-4 

1529.59 

1500-0 

16292.9 

10561.1 

12.5 

1531.27 

1600.0 

13865.9 

11774.6 

9.6 

1532.94 

1700-0 

10908 . 2 

13253.4 

6.7 

1534-62 

1800.0 

6755.3 

15329.9 

3-0 

1536.29 

I 

1 

J 


2.5.2  Plots 

All  graphical  outputs  are  plotted  on  the  Tektronix  Terminal  and  have 
a  similar  format.  The  results  are  plotted  on  a  5x5  grid  and  a 
scale  factor  for  each  parameter  is  printed  on  the  display;  this  scale 
factor  is  the  number  of  metres  equivalent  to  one  division  of  the 
display  (i.e.,  full  scale  is  five  times  the  scale  factor).  At  the 
end  of  each  display,  a  copy  is  produced  on  the  Hard  Copy  Unit. 

The  origin  of  the  plot  is  the  upper  left  hand  corner;  the  vertical 


1  yr*f  ••  .J5  <**»l 


■  J%Ur*>~-'Z  fc'-  **«*,  ‘  • 


parameter  is  depth  and  the  other  relevant  parameter  is  plotted 
horizontally . 

a.  Distance,  extent  and  naxinun  thickness:  one  plot  is 
produced  for  each  parameter.  Examples  are  given  in  Fig.  5- 

b.  Shape:  One  plot  is  produced.  If  sore  than  one  source 
all  outputs  are  plotted  on  the  sane  graph 

An  cxaople  is  given  in  Fig.  6.  \ 

| 
j 

1 

I 

! 

I 
is 

load  the  data  tape  in  the  tape  reader  and 
enter  the  program  (under  the  Real -Tine  Systeu,  this  is  effected  by  j 

typing  “OX,  SHADZO “  on  the  Operator  Console).  If  the  last  point  ] 

* 

i 

on  the  tenperature  or  sound  speed  profile  is  not  deep  enough  to 
allow  the  critical  depth  to  be  calculated,  the  conputer  prints: 

THE  I  AST  POINT  OF  THE  B.T.  IS  ABOVE  THE  CRITICAL  DEPTH 

and  the  prog  ran  ends.  In  this  case,  the  data  tape  nust  be  retyped. 

Once  a  correct  data  tape  has  been  read  and  the  sound  speed  profile 
has  been  calculated  and  listed  as  shown  in  Table  2,  the  conputer 
will  print: 

FOR  SHADOW  ZONE  SHAPE  TYPE  1  OTHERWISE  TYPE  -1  * 

2.6.1  Distance,  Extent  and  Maximum  Thickness: 

If  -1  is  replied  to  the  question  about  shadow  zone  shape,  the 
computer  outputs: 

FOR  GRAPHIC  OUTPUT  TYPE  1  OTHERWISE  TYPE  -1 


depth  is  requested, 
to  the  sane  scale. 

2.6  Operation 
To  run  the  program. 


All  questions  of  this  type  should  be  answered  by  either  -1  or  1,  followed 
by  carriage  return,  line  feed. 
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and  after  the  response  to  this  request,  the  cosputer  outputs: 

FOR  MULTILAYER  TYPE  1,  FOR  2  LAYER  TYPE  -1 

When  the  type  of  sodel  has  been  selected,  if  no  graphic  output 
was  requested,  then  the  coaputer  prints: 

TYPE  IS  S. DEPTH:  START,  STEP,  END 

The  source  depths  for  which  output  is  required  should  then  be  typed 
in,  giving  the  ninisun  depth  first,  and  with  values  separated  with 
comaas.  So  output  Is  produced  for  source  depths  greater  than  the 
critical  depth.  When  the  source  depth  has  been  input,  the  results 
are  calculated  and  printed. 

If  graphic  output  is  selected,  then  the  three  graphs  are  displayed 
and  a-Tter  the  last  output,  the  joystick  is  enabled  and  the 
computer  prints: 

TO  EXPAND  PLOT  POSITION  CURSOR  AT  MAXIMUM  SOURCE  DEPTH  REQUIRED 
AND  TYPE  1  OTHERWISE  TYPE  0 

If  an  expansion  is  required,  then  the  horizontal  line  of  the  cursor 
should  be  placed  at  the  maximum  source  depth  for  the  new  plots 
and  1  should  be  typed.  A  new  set  of  li splays  is  then  produced, 
and  the  expansion  option  is  repeated.  (It  Is  not  possible  at  this 
stage  to  increase  the  naxiaun  source  depth  of  the  display  -  If 
this  is  required,  the  progran  must  be  continued  and  reprocessing 
of  the  same  data  request  ed.  See  below) .  If  expansion  of  the  plots 
is  not  required,  type  0. 

2.6.2  Shadow  Zone  Shape 

If  the  shape  is  requested,  then  the  source  depths  are  requested 
as  for  the  printer  output  above.  The  source  depths  must  be  below 
the  depth  of  minimum  sound  speed;  if  the  starting  depth  is  less 
than  this  value,  then  it  is  automatically  adjusted  to  this  value. 
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If  the  source  depth  exceeds  the  critical  depth,  the  calculation 
stops.  If  only  one  source  depth  is  required,  the  second  and  third 
parameters  should  hi  zero.  The  computer  then  calculates  and 
displays  the  shape, 

2.6.3  Termination 

After  the  relevant  output  las  finished,  the  computer  prints: 

TYPE  0  TO  REPROCESS,  1  TO  FROCEL  '  MEW  DATA  Ok  2  TO  STOP 

If  the  response  is  0,  the  program  returns  to  the  point  at  which  it 
asks  if  the  shadow  zone  shape  is  required.  If  the  response  is  1,  a 
new  data  tape  is  read  and  processing  restarts.  If  the  response 
is  2,  the  program  ends. 


it*  ~X*r«  **. 


nr  Ft  u 

CRITICAL  DEPTH  *  1562,15  METRES 


^  DCTTH  SCALE  *  20  M/tSIT  RAKE  SCALE  *  10000  M/CSIT 


IIP 

jggg 

WM 

\y 

_ 

SOURCE  DEPTH  =  100.0  M  TO1700-0  M  IS  STEPS  OF  200.0  M  XAHGE 

j 
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FIG.  6  EXAMPLE  OF  SHADOW  ZONE  SHAPE  PLOT 
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APPENDIX  A 


DETAILED  PROGRAM  INFORMATION 


A.l  Flow  Chart 


A. 2  Explanatory  Diagrams  and  Program  Listing 


Note ;  The  word  "velocity"  used  in  this  Appendix  corresponds 
to  the  words  "sound  speed"  used  in  the  Main  Text. 
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PAS £  3*31 


FT3.L 


T-siS  *£S  DEVELOPED  AT  - 


S*Cl*4T  as*  research  centre 
VlALE  S an  SARTOLG'CO  -  A## 
l*52fe»  LA  SPEZIA 
ITalt. 

Snji22(3»7ii 

J.P.  27-14-71 

OlSS&SIGn  d(44)tT(^)  »CH«)  *1  ASCI  (19)  «INAH£(6) 
DIMENSION  «A (1331  «£X ( lv2) 

THIS  ?0U'j3AH  PREDICTS  T-mE  RANGE  AND  EXTENSION 
OF  THE  Sh*03«  ZOfcE  FOR  VARIOUS  SOURCE  DEPTHS. 

READ  DATA  TAPE  a?iD  CONVERT  TEMPER  AT  UkES  TO 
VELOCITIES  IF  NECESSARY. 


PHI=42. 

SAL=38.6 
13  00  23  1=1.6 
23  iWkHECl) =233434 

READ (5. 33)  (Ir|Arf£( I)  .1=1.61 
30  FORMAT (BA1) 

REeD(5«*)lDAY.H0NTH.lYEAR*lH0JR.>*lNUT 

READ (5. *1  MODE 

fc=I333(:<00£) 

00  35  1=1. N 
35  READ(5.»)DfI>fV(I) 

IF (NODE) 43.53 

43  CALL  3TrtA*l(V.0*>J*SAL.?HI) 

50  CALL  FIHTP(V.N.ITURN) 

0TORM=C ( 1TURN) 

VTu.RV=V<ITUkN> 

OUTPUT  HEADING  AND  VELOCITY  PROFILE 
CALL  EXE C(3*11063.-l) 

WRITE (6.63) (IHAME(I) *1=1.6) *IOAY. MONTH  *  I YEAR . IHOUR • MlNUT 
63  FORMAT (24X.".SHAD0.v  ZONE  P*EDICTI0NS,,»//.7X."XBT  ".6A1, 

1  12* *UD-  TE".  13.2  12)  *12X*,,TIME,,.I3*,,-"*I2.//. 

2  5X  *  "DEPThS"  • 1 0X . "  VELOCI T I  ES"/5X  .  "METRES”  *  10X  » 

3  "METrtES/5 rC"*/) 

WRITE (6* 70) (0(1) *V (I) *  1=1 *M) 

70  FOrf'-SAT  (A»F10.2»rl/.l) 

FIND  CRITICAL  DEPTH  ETC. 

IF(V(»)-V0)d0.100 
80  WRI T£(l«9o) 

90  FORMAT ("ThE  LAST  POINT  OF  THE  B.T.  IS  ABOVE  THE  "* 

1  "CRITICAL  DEPTn") 
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PAGE  m**Z2 


C 

C 

C 


C 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


STS? 

V2=*;«lj 

4'=?. 

IF=1 

Co  12*  l=2:.i»i 
iri^rMi)J!Wf!24 
11t£  ~*='-»iS 

|?j  eS*sTS.Vj£ 

C^LL  ”1  ®L (’/ *  ICftlT *N«lTUiM) 

VCHIT**«S 

5RITE  Jfc  *  1341  OCfrIT 

13tf  F€-a^TC/^12*»“CSITiC*L  OEPTrt  =",Fd.2.  »  METRES"/) 
tfEuuE* T  TYPE  jr  CALCULATION  AND  OUTPUT 

1^4  -PITEil? 

l**r  ~  -  ,m*T  ji  2uvE  SHAPE  TYPE  1  OTHERWISE  TYPE  -1 

*££351**1 15H4P 
IF{t^5Pli54v72» 

15r 

155  ("rO.'<  OHAPmC  OUTPUT  TYPE  1  OTHERWISE  TYPE  -1  «-"l 

R£»3I1.*1IGkAP 
SPITE  Ciil-52?i 

163  FvP-ssrC'FoP  hjlTIL«Y£P  TYPE  1,  FOR  2  LAYER  TYPE  -1  «-»•) 
K-£0!i»5)LiTt(i 
IF(IS3APJ  i7»i*35a? 

PslNTtK  OUTPUT 
17-4  iPiTEd^ls^l 

lei:  F~i2*&T(“:YP£  In  S.i/EPTri:  START, ST£P*£nD") 

RE-  J  ( i  •  )  sO  i  »  So 2  *  Su3 

OUTPUT  rtciJlMu  ETC. 

CiLL  £/£C(3?ili)6S*-l) 

IF  (LAYER)  ISi’,213 
199  *‘S-'IT£(5«2a*> 

2a-5  rC-i-sAKcSA.-'HO  LAY£«  MODEL"/) 

GO  TC  23(5 
21  i  !* .91 7E  {-■:  ,22-5) 

25 f yP  SA 7  {  2.yx  .  » m*Jl  7 Ila YEk  MOuEL"/) 

23/.*  .PI  ic.(c,2-*'j) 

24r  FoR-iATC*  SOURCE  D£rTH‘,4X“£XTENSIGf*l',7X“R4NGE"6X"THlCKNESS" 

1  5X»“S.  v glOCI T r •*/ . 4X , 11  ME TRES" ♦  t*A , ‘'METRES" , 9X , "METRES" , 

2  6X. "METRES",  l«M,"r1/S£C") 

CALCULATE  SOURCE  DEPTH  ETC. 

Nx= ( 333-50 1>/S02*1. 

oc  y-/-  l=i  ,\x 

SlEj=S'  i ‘•FljA r ( I-i i  -su2 


C 

c 

c 


c 

c 

c 


IF  OJrsS&l  oif«344 
IFCSSiy-liCsingaSjj^ 

26«*  I*"  i  SOZP-UTite  i)  27*124* 

27cS  CJLL  rr<AL«a.4iSu£f-,*tL*I!Jti»jlT^;i) 

5f  L~  /  *  JJL  1  »“*3*  «0L«9EP  *  I9.M.  ITtMfil 

GO  13  29» 

29i?  Si!’!’ 

SfKL  rI*AHv.O,»Ei  *3tP,Iv*I&,ITU«iij> 

silt 

T*0  L*Y£s?  aOfla 

3M  ^Li0r32*V*‘VC!<I  r*"“»I*^w.ori^.vrU«.«t^,E,T.THia. 

aULTlLfcYER  :-500£L 

31a?»ic<S"''**V£L*,C"Ir*0“,T’S3eP-IC«T-IW’V,0,R»„GE,EXT^ 

PrtIf4T  «£S-JLTS 

S  rf£If i6-?301 59EP  •  t  a  t  *  SA*3E  *  thick  *  yel 
334  FO^HATfr laf.l,3Fl4.1*FlS.21 
342  CONTINUE 
GO  TO  15z2 


C 

C  »5'H}-553.355iS*v««45SS,j^jW55jf4s 

C 
C 
C 


GS4PHIC  OUTPUT 


35SJ  SOM  A  A=«iCfc  I  T*iJ  .95 

SOhINs  (5U/J4A-0*;)  /i*4.  *00 

FlftD  vlAAl***;M  VALUES 

36tf  IF (LAYtH) 374*333 

37*  SA.LL  J,u*v'  ’“Crf*r*£'C«IT*S0KAX,03,0rUrtM*tfTU«M.»AMAif  c*  TJ1 
GOLrOT^y"*^CKlT’UCKIT’SOHiN,0a*OTUKN*VTU,iN*RA*EXMAxlTHHAX) 
3«5J  IF  <SD9*a-DTUR*I)  394*4*3 
3-*vi  call  Fl7AL(0*V*5D'i/,X*y£L*I0EPf  jtu*iV*I) 

CALL  r I V«L(V*0* V£Li JtP* I9tW* ITURNI 
GO  TO  <tIv 

Call  ^IVAL(0*/,5C.MAA,VtLtIDEP,N*ITUKN) 

CALL  r  1YAL(V *0* VcL*0£P*  19*  14* I TURN) 

CALL  F  i  t/AL  ( D  *  v  »  OcP  *  VfcL  *  1 7  *  I  TURN  *  1 ) 

41.,  “KL  n^UV.O.VEL.^.I^i/.ITUm. 

*  v  9 ij  f  *£,4; 

CALL  FfVAL  (£*'/*  SOXi.  9  *  VEL*  I  l»£P*  I  TURN*  1 ) 

CALL^r  I VAL  ( «t>*  vEL»0£P»Iv*N*l  TURN) 

1  “L%^K«!EXM«),iCKlT’UCRlI’SDM1H,ICBIT’I0EP* 
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CALL  -OARA  (624*644) 

C*U  ALPHA 

iFl J— 2) 43ar*464*49tf 

432  af^iT£  ( 16«44<£}  ISCTm 

FC^«Tt“T!jICK<€SS  SCALE  =*,I5t»  H/4MIT") 

CALL  DaS!K  1912, 151 
CALL  ALPHA 
•RITE (16*454) 

45e  FORMAT  (** rnICK’JESS«-Ml 
GO  TO  523 

463  5~ir£<  16*470) I5CRA 
47*i  FORMAT  SCALE 

CALL  0ASIUO12.15I 
CALL  ALPHA 
WRITE (16**#9) 

43.3  FCwWT (,*fiA-JSE<-I,l 
GO  TO  523 

493  VPI  T£  ( 1 6.SZ/5)  ISCEX 

533  FORMAT ("EaTENSION  SCALE  =M*I5*M  h/UNIT") 

CALL  0***1912, 15} 

CALL  ALPHm 
WRITE (16*510) 

510  FORMAT ("EXTENSION*") 

529  STE?=  (S0HAX-D3) /YZti, 

C 

C  OUTPUT  RESULTS 

C 

00  553  1=1. 1-7.) 

S  D  E?=D  ?-  *r  lO  -T  ( I )  -STEP 
IF (J-2) 525*533 
525  IF(LATEk) 533*540 

539  CALL  T!>=O(vv'*VC*iT*«JCJ<IT*SOEP*D0*DTuRN*VTURN* 

1  *A(I).EX(l)*Tn} 

GO  TO  5taa 

549  IF (SOEP-OTURH) 550*550 

550  CALL  rIv6L(U*V*SOEP.VEL*IOEP*ITURN,l) 

CALL  r  1  y al  ( v  *  D  ♦  VEL  *  0£?  *  1 9  ♦  N*  I  TURN ) 

GO  TO  57-3 

569  CALL  FI VAL (9*  V*SD£P* V£L*ID£P*n*ITURN) 

CALL  r I VAL  < V ,0 ♦ VEL *DEP* 1 9 *  10* I TURN) 

CALL  FIVAL(D*V*DEP* v£L*I7* I TURN* 1 ) 

C A L L  F I V  AL ( V ♦ D  *  V  EL  *  OEP  *  I d » 1 7 . I TURN ) 

57*  CALL  MULTI (V0*V£L,VCKIT*OCRir*SOEP,lCRIT*IuEP,V,O* 
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PASS  ^7-^5- 


1  ft*fIJ-£*CII) 

CALCULATE  CS-u*Oi:iAT£S  FC*  sSISPLAf 


53s*  Ir  <  J— 59it3t<«6l  J 
59*  U=( W/rU>6H liCIrtl J*<J@*.*23. 

60  TO  62» 

6*2  I^J^ai/FtOirdSC^D^iW.^. 

GO  TO  62® 

614  IX=I£a  « 1 1 /r  OIAT C  ISCEa)  I *23** *23. 

624  IYs»31.-CS2£>/FLS«I I15CSOI 1*12®. 

lr  f  i— 2>63«<S*r»43 
63»  C* LL  3A»<KIIX*Ifl 
GO  TO  6Ss» 

64«  CSLL  5$iG--sSlX«m 
654  CGmTIViE 
CALL  HARDC 
562  C*^TI»j£ 

Call.  =j*tsK*?.<?«331I 
CALL  ALPHA 

REQUEST  IF  cXpMiSIO*  IS  REQUIRED 

*»3IT£(16«S74J 

67»  fOr'RATCTO  EA--4-iii  r-LUf  PtJSlTIGH  CURSOR  AT  MAXIMUM  SOURCE  "* 

1  "CEPTw  H£*juI«tcO  A69  TYPE  l'**/*3i*X,‘' OTHERWISE  TYPE  3»«l 
CALL  CuRSI (ICHAR.Ia.iYJ  C 

I r  ( IC-iAR-491  71»«6-3» 

CALCULAie  k£#  .-SaaIHuM  SOURCE  DEPTH 

63s?  S0MAX={fc3i. -FLOAT (J.YI  l/120.*i'L0AT(ISCSD) 

IF  (5D'JAA-uCr<I  /-5**69w 

694  S0MAX=uCf<lT*»*v5 
7-40  GO  TO  36*.- 
714  CALL  z^Abc 
CALL  HOME 
GO  TO  1504 

vv*>v^vvvvC'GviJvvv,i/vvvvvv^rif,5  lUfr-Jr^KHrifr-sJ^fr 

CalJuLajIO;*  A.-jJ  oISPLAY  Or  SM4U0^  ZONE  SHaPE 

7 20  r'PITE  ( 1  *  Is?) 

REA  j  { 1  •*)  :>01  •  502*5l>3 
IF  (S:Jl— DTiJ''  •»—  i  ■ )  7 5«>«  744 
734  SDl=OTUH.'i*i. 

744  IF  (502)  7b^  •  75-?  ♦  764 
75.'  ^-<=1 

50  TO  7ab 

Tty  IF  (SOi— bUi )  77v)*  774*  1*i 6 
77-3  Au  =  l 

Gw  fu  Tob 

756  H^  =  t5-j3-iJ3  J/SL2+1. 
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puce 


FlfcD  MAXIMUM  VALUES 
7-35  S5£?=S©  1 

CALL  FlVAL(U«V*SO£P«VSOEP*ISOPT*N*ITlMfi» 

CALL  FI VAL ( V*0* V50£?«DMAX* I * I»« ITUR*» 

CALL  rlvALf£r«V*$HAA«Vi»Il*ITu4NilI 
CALL  F I  VAL  ( V*3*  tfhulM*  1*11  «I  TURN) 

DC  79i> 

799  C(I)=(«j(l)-u(i-l)?*(V(i>*V(I-l>) 

Z=0«^)/l«i  * 

CALL  U*d««Stf cP*  VSUcP*  ISDPT*i5»«X*  19*  ITURH*N* 

I  C*Q.V*uISri*uI5T2*kA0) 

CALCULATE  SCALE  FACTORS 

ISCTh=ISC*L  CsJHArAl 
iSCRA=ISC*L(OIST2I 

OUTPUT  bRAT  IC'JLc 

CALL  3P*T(  iiTAME*  ISCTrf»OC«JIT*LAY£R* ISHAP) 

CALL  DaRK (633*544) 

CALL  ALPHA 
WRITE (16**72)  ISCRa 
CALL  DARK («.'3) 

CALL  alpha 
IF  (VX-2) 892*912 
89i?  SSITEdi^^alSOl 

900  FOR *AT ( "SOURCE  DEPTH  ="*Fo.l*"  METERS  *«) 

GO  TO  930 

910  H*ITE(lb*923)Sul*SG3*SD2 

92£  FORMAT ("SOURCE  DEPTH  =”*Fo.l*"  M  T0"*F6.1* 

1  “  *=  I*  STEPS  UF"*F6.1*«  .M«-") 

933  CALL  DARK (912? IS) 

CALL  ALPHA 
WRITE (16*480) 

DO  11  oz*  J=1*nX 

calculate  RESULTS 

SD£?=SOi *F LOaT ( J-1 ) *SD2 

CALL  FIVAL(D*7*3D£P*VS0EP*I5DPT*N*ITUrtN) 

CALL  F L vAL ( v  O*  VSutP*uKAX *11*10*1 TURN) 

CALL  F  i  vAL  (D*  V*u<'jaa*V1  *  1 1  *  I  TURN*  1 ) 

CALL  F[vAL  (V*0*V1*'>'AX*I2*I1{.ITURN) 

STE3=(G MAX-03) /ISM. 

ZOLDlss*. 

Z0LD2='J. 

RNEtfl  =  .«. 

RN£*2=1.E19 
RAD 1=3. 

RAD 2=3. 

00  1 351-’  l  =  l*l*-3w 
Z=FL0AT(I)*6TEP 

CALL  0rTPS(Z*0D*SDtP*VS0EP*ISDPT*UrfAX*l3*ITURN*N* 
1  C*r>*V*0ISTi*0I5T2»K4D) 
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C 

C 


F im  •HlZti  POINTS  TO  RETAIL 

IFf*'A0iT9=»2*'i>4s}*9539  I  ** 

60  TO  962  ** 

95^  E=(Z-2lLU1)/RA;>1  « 

RA*g5^aE*l*34itf  «St{4T  (E*  (2.-E)  I 
96J  IF  (3IS  «9di 

972  3a(I)=F;A?J5£ 

GO  to  -93 

932  s?-:£-l=Vr2tl 

R4CI)=0I5T1 
Z0L31=Z 
RAs)  1=3*0 

99fc  IFCRAD2J  l^tlsa^tl/lB 
Hi-sSE^-.^EirZ 
GO  TO 

1919  E={Z-Z0L02)/kAD2 

R4';5E=^&P2-R4Q2*5-j-cT{£*(2.-£) ) 

I32m  IF  { R  A'*GE-0 1 S  T  2  )  123 3*  i340 
133«  E/.(I)=^AfJGE 
GO  TO  1X52 
1342  R>J2£2=SISr2 
EACD-0IST2 
Z0L32*Z 
RA02=-»  Ai) 
l/f5i?  CO^TI'juE 

DISPLAY  RESULTS 

00  U’S-‘  1=1  *  l'.ivJ 
Z=FLOi Til) *bTEP 

I/.=  (f=-  ( I )  /FlOAT  ( I  SCR  A) )  *299. *23 • 

IY=s31.-{Z/FL0aT( 1SCTH) )*122. 

IF  C 1-2)  126-3*  l'.l72 
1263  CALL  DARK (IX* I Y) 

GO  TO  1333 

1373  Call  SPIGrt  ( I X*  I Y)- 
1-330  CO.-jTI'iUE 

00  139/  11=1* 132 
1=151-11 
Z=FtO*7 (I) -STEP 

IX=  ISA  < I ) /FLO*  T (ISCRA) ) *232. +23. 
lY=331.-(Z/.-L0^T{IbCrH)  )*123. 

I- 392  CALL  S.-'IG-jUX.IY) 

II - 33  COwTI'juE 

CALt-  -iARDC 
CALL  E^ASE 
CALI  ALP^A 
CALL  HOME 
GO  TO  1523 

-2  -5*  •»«***«•  «■»»•»**■»**•»*# 

REQUEST  NEXT  OPERATION 


page 


c 

15^i?  •’inlTz.litlolc) 

151*  FO-<SATC,,Trt»£  *  TO  PE^-cCCESS*  1  TO  PROCESS  NE»  DATA 
1  "0^  2  T3  STOP") 

R£aD(1«*> i 

IF  ( I-i  )  !>/  •  l  152* 

152*  STOP 
END 
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J.P.  22-U-71 

d- i:>  ,!f1C  FI“(ST  S£T  0F  SUBROUTINES  FOP  THE 
•  Srt*20 

SLiu-<0JT1nc.  '-OLTI  ( V/|«  VEL*  vCkIT»DCRIT,SUEP»ICRIT>IDEP» 

1  V*L',h«NGc»EaT) 

0I-M£  iSU*  v  ( 1 )  *  j  ( i )  ,  T  (42? ) 

,  b,^r0jriNt  Calculates  the  range  and  extension  of 

T  .c.  5  XV..  Eu-it  F 0*  1  HE  MULTILAYER  MODEL • 

OCR  1=0(10  IT) 

VCriI=V  ( ICR  I  D 

oo  io  i=''*icRir 

VEL 1  =  V (1DEP) 


calculate  extension 


If  (DEPl-l'CPI)  23,3a 
20  I  =  I0£P+1 

of  p2=d  ( i ) 

VEl2  =  v ( 1) 

*  *  vkJ-vEL)  )  +  SQRT ( ( V0+VEL2)  *  ( V0-VEL2) ) 

exTa(JtPE-S0aP)*<vtL*VEL2)/S 

S=S'JRT  i  (  VM  +  vCkI )  MvA-V  CRI)  ) 
tXT=EXT+  (L’CRI  T-UCRI  )  (  VCR  I T+-VCRI )  /S 
GO  TO  4 A 

ll  i m«^;mMr<'T*'a,/sa,T,w,V£U,,vj-'a" 

50  11=  IOER+2 

DO  SA  t=Ii,ICRIT 

60  ExT=tXT+(D(I)-U(I-i) )*(v(I)+V(l-l))/T(I) 

7 w  EXT=fXT*2, 

CAu.^ULmTE  ^ANGE 

S=S jRT ( (  VtLl ) * ( vA-vELl ) ) ♦SORT ( ( V A* VEL) * ( VW-VEL) ) 

RA  <5-  =  {  Su'cR-uEPi  )  *  ( VtL+VELl ) /S 
00  5  a  1=2,1 u£P 

R  A.m3E=-.Ai'JoE>  (J(I)-0(I-l))*(v(I)*V(I-l))/T(I) 

R  c.  T  J  ">  i 
END 
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SUSROuTlNt  T«O(V'0.VC'ViT*DCRIT*SDEP*O0*DTURN* 

1  VTuR«»  *>aNG£  *EaT*  THICK) 

THIS  SUBROUTINE  CALCULATES  ThE  RANGE*  EXTENSION  ANO 
ThICaVcSS  OF  TnE  SHA'JOw  ZONE  FOR  THE  2  LAYER  MODEL. 

IFOT  jrfN— Si'Er’) 

10  V£L=tf  T’jrth*  ( vCRI 1  -v  I  jkn)  *  ( SuEP-DTUKN)  /  (DCRI T-OTURN) 
S=S0RT(  <V0*V£L)*(V0-V£L) ) 

EXT=2.* (DCPI T-SOEP) * ( VCRI T* VEL) /S 
S=S*S3~T ( (V0*vTurn)*(V0-VTURN) ) 

RAN3E= ( SOEP-OTURN) * ( V  TURN+VEL ) /S 

RAn3E=PANGE-*-OTUKN*(Vki*VTURN)/SQRT<  (V0+VTURN) *<V0-VTURN) ) 
THICK=03* (OTURN-D0) * ( V0-VEL) / ( V0-VTURN) 

GC  TO  30 

20  V£L -VTvjRfc*  ( V0-v  Turn)  *  (SOEP-OTURN)  /  (O0-OTURN) 

S=SQRT ( ( Vt *  V TurN) * ( Vw-VTURN)  ) 

RAN3E=SDEP*( V0*VEL) /S 
ExT=2.*(bCRIT-DTURN)*(vCRIT*VTURN)/S 
S=S*S0PT  (  ( V0+VEL) * ( V0-VEL) ) 

EXT=EXT-2.* (SOEP-uTURN) * ( VEL+VTURN) /S 
ThIC<=OTURNMDCkIT-DTUKN)*(VEL-VTURN)/(V0-VTURN) 

3>J  RETJRN 
End 


E  0003 


SUBROUTINE  3TRAN(V*DfN*S*PHl) 

DIMENSION  V  < I ) .0(1) 

THIS  SUBROUTINE  CALCULATES  THE  SOUND  VELOCITIES  V(I) 
AT  DEPTHS  D ( I )  FROM  THE  TEMPERATURE  VALUES  INITIALLY 
STORED  IN  AncKAr  V.  S  IS  THE  SALINITY  AND  PHI  THE 
LO  -1301  ruuE 

DO  1 0  I  =  1 » N 

V  ( I )  =l'.v3.*3.«(V(I)-l0.)-J.006*(V(I)-10.)*(V(I)-I0.) 

1  -K.04'MV<I)-lB.)ft(V<I)-18.)*1.2«(S-35.)-0.0l* 

2  (V{I)-lrt.)*(b-JS.)*D(I)/M.*l.£-7*D(I)»0<I> 

3  ♦ 2 . £- 1 /*U (i)*J(I)*(V(I)-l8.)*(V(I)-18.)*l .E-4» 

4  0(1)  - -Hl/Aii.+c,  oE-4»  V  ( I )  *  ( V  ( I )  -5. )  *  ( V  ( I )  -25.  i 
1">  CONTINUE 

RETJP\J 
End 
c 
c 
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Sue ROJ Tint  F I N T  P  ( v  *  M » I  Tuk.'J ) 

DIMENSION  V  C 1 > 

THIS  S'JaROuTlNE  FINOS  TrtE  i*5 1 N  £  MUM  ELEMENT  OF 
ARRAY  v  AND  PLACES  ITS  SUFFIX  IN  LOCATION  ITURN 

VMI N=V ( 1 ) 

oo  i,f>  i=e*n 

I F  ( V  { I )  - v P I •  i )  1  u » 2  u 
10  V*IN=Vd) 

IM£N=I 
20  CO-iTINuE 
ITUHN=lMlN 
RETURN 
END 


AGE  0005 


SUBROUTINE  FI VAL ( A»3 » VALUE *  RESLT * ISSCR tNl *N2) 

Or-tE  lSION  A(l)  ,3(1) 

THIS  SUBROUTINE  FINDS  TrtE  VALUE  IN  ARRAY  B  WHICH 
CORRESPONDES  ru  "VALUE"  IN  ARRAY  a  and  PLACES  IT  IN 
RESLT.  Nl  ANO  N2  ARE  THE  LIMITS  OF  THE  ARRAY  SUBSCRIPTS 
SO  ARRANGED  THAT  N2  IS  THE  SUBSCRIPT  OF  THE  LARGEST 

value  of  a. 

N=I A3S (NI-N2) 

IF (N1-n2) 10*20 
10  INCR=1 
GO  TO  30 
20  INCR=-1 
30  OO  50  Il=l*N 
I =N 1  +  INCk*  1 1 
J=I-INCR 

IF ( VALUE- a ( I ) ) 5U. *0 
40  DEL= ( VaLUE-A ( I ) )  ‘(A(J)-A(I)) 

RESLT=H (I; +DEL*(3(J)-B(I) ) 

ISSC»=I 
GO  TO  60 
50  continue 
60  RETURN 
END 


onor.  ooonoo 


-  ’  -1 


FTm.L 

J.-.  27-iWl 

T- 1  b  IS  T 8c  5£.CJ:jO  SET  0r  S JEhO'JT i?-.ES  FOk  THE 
S-.i'iQ  " 

SlB-'JJTI-.E  GR  *7  ( I..*»  -ie.«  I  SCdD»uCP  I T ♦LAYER*  ISHAP)  : 

0!  'E  5  I  -•  [  !-a  (i) 

THIS  SUdHOUf I:«i  pkuO'JCES  ThE  HASIC  GRATICULE  FOR 
GRAPHIC  DISPLAYS. 

CfiL.  E-hSc 
CALu  HO’-it 
WRITE ( 1 6» i ") 
lei  FOR  T  ( /  ) 

Ca;__  DoUtlL 

IF  ( I  5 n -» r )  —  » 2 J 
2A  WRITEdb.Jd 
3  vi  F0H'1dT  (2^  A  •"*•") 

40  VnR  I  TE  ( 1 6  » 5  £ )  ( 1 1  j  4M£ ( 1 ) *  I =  1 »  6 ) 

So  FOn^AT  l  12a.»'X  d  T  "S6(X«A1)  >"<•") 

IF  (LAYER)  D-  *3-C»l;i.J 
60  '.*/ RITE (15.  7-i ) 

70  FOK'IAT  (  22a  ♦  ‘T  ,v  0  L  A  Y  E  R") 

GO  TO  1 20 
80  WRITE (16*4/) 

9w  FORMAT  (A) 

GO  TU  120 
15)0  >^I  TEd6.il  d) 

lio  Fjrviat  (EEa  ."  -I  j  L  T  I  L  A  Y  E  R") 

12')  CALL  AL-- 

W*ITEdo»13o) JCRIT 

130  FOR  HAT  ( /  »  27  a  ♦  "Cr  I  T I  CAL  DEPTH  =".Fd.2»"  METkEo1') 

IF ( ISJ -P) 1-if 160 
140  WHITE ( lb»lSo) 

155’  FORMAT  ( oX . "SOURCE  «■") 

GO  TO  160 
160  WRITE ( 16.170) 

170  FORMAT  (  13a  t  »«■••) 

18/  wITEdb.i-  )  ISCSO 

190  FORMAT  ("uEPTh  SCALE  =  "*I5»"  M/UNIT") 

CALL  DARK  (23»->31 ) 

CALL  3- 10- (23.31) 

CALL  3RlOn(1023.31) 

CALL  3RIbn ( 1023.631 ) 

CALL  9^1  On (23.631) 

DU  210  I  =  223  •  823  «  2‘,15' 

CALL  DAkK ( I . 31 ) 

CALL  3R Ibn ( 1 . 63 ) 

DO  2 t'i  J=c3*563»16 
23/  CALL  ?OiwT(i»J) 

CALL  DAkk(I*S94) 

21/  CALL  9>IGn(i.o31) 

DO  23/  I=i5i »51 1 . 123 
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PAOE  i? 


CALL  3aKiC(23«I> 

CALL  3-/lGn  (55.  I ) 

00  22/-  J=a7. 991.1b 
22  ’  CAL'-  ( j.  I ) 

CALL  3akM991.I) 

230  CALi-  lt>23»  I ) 

CALL  0ahK(7.o31) 

Call.  *Lr*rt*» 

9SI FE ( l6«24w) 

24y  FQ-H4T  <••«?••.//) 

IF  ( IShAP) 252.270 
250:  i»PIT£  ( l&*2*>c) 

26vi  FOk^at (••s**/,*o-,/"U'*/MK,,/,,C,,/,,f*V) 
270  vRl FE (16*280) 

28M  ,F0H*«T  ('»U"/*‘£•'/,,P•,/*•  ^"/,,rt•,) 
RETURN. 

END 


AGE 


! 
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Eu  iCriO't  1SCAL  (kAhaX) 

01  'lE''JS  IUN  ISC  (3)  »/HAX  (3) 

i  this  function  iakes  r.-it  mexf  value  above  ramax  from 
1  hE  SE1' ic.b  2tt  1  V7**i\i *  5*10**N  WHERE  N=1.2*3 . 

ISC( I)=i 
ISC (2) =2 
IbC ( 3 )  =5 
XM«X ( 1 ) =5* 

XMAX (2) =14. 

X  w  A  X  (  3  )  =  d  z  » 

1  «d’  00  3 i  1  =  1.3 

IF  (  X'^AX  ( I )  -RAMAX)  3x3.2-!) 

20  ISCAL=ISCU) 

GO  TO  50 
30  C0nTI\jE 
00  '■*'!  1  =  1.3 
.  lSC<I)=l0*IbC(I) 

40  XVIhA  (I  )  =10.*X  AAA(  I  ) 

GO  TO  lw 
50  RETURN 
END 
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ooo  non  o  n  r>  non  nnnnn 


?*'. j£ 


So'3; )« I“:i  j-iT.-i(ZoLi;P»^SDci,*I5D',f«D,l4X»IiJ*ITUSS*N* 

l  c.c  *v»«ji5n  »oi:>rz«RAO) 

DI-Z'.SIO-.  C(l)  *001  *v(l) 

T-ilS  5*J3-<0:jT I-.E  CALCULATES  Th£  DISTANCES  TO  THE 
VC-Tc  -£5  OF  KAYS  -niCi  OP  T  £  A  AT  a  OEPTH  Z  AND  ALSO 
THE  P4ZIUS  uF  Tn£  nAT-  AT  THESE  POINTS 

ps  :  u. 3  >=:>*;«  M  c**3j*<  a-S)  > 


Flu  vELOCiTY  mT  JEPfn  Z  ETC. 


DS  33  T=OU  SOP  T) 

VS-->->T=v  ( iaOr'T) 

CALL  FlVAL(U*V*Z»:>A*iZPT.i  TitoN* 10) 
IZPr^iZHf.i 

DZPr=0<IZPT> 

VZPT=V(iZPT) 

SX=Si'*:£LLS  CO.iSTanT  FOP  THE  PAY 

CALL  "I  --L(  /0.o'.^_Z.ILZ-T*N.ITO.P-J) 
ULZ*T=.>(IlZPT) 

VLZJT  =  v(I!.ZPT) 


fi lj  fnc  oistmwlc.  to  the  fip^t  turning  point 


oisti= 

IF  ( .■  S  „  r  \ )  3  »  i 

3  01  ST  1=  { Su£C—uS0P D  *  *  VSDEP+VSDPT)  /  IP 50  (SK *  VSOPT)  ♦RSQ(SK*VSDtP)  ) 
7  IF(IS:)?r-I2?T)l4«3»'l4 

00  2)  I=J»IbOH,f 
IFtSfWd-I)  )  17*13 
13  IF  ( SK-./  ( I )  )  17*20 
17  L1  I  r  i  .  —  • 

01 5T2=  1 . £  1/ 

GO  TO  1A0 

20  0I3Tl=0I5Tl+C(I)/(r<5U(SK*</(I-l))-*-KSQ(SK*V(I)) ) 

30  IF(l.Z-4-SK*\/ZPT)4o,45 

4,  : 1 371=015  n*  (OZPT-Z^IVZ^T  — ...  /r3u(SK*  VZPT) 

FI:|u  T.-iE  QISIaimCE  TO  Tn£  St-COND  TURNING  POINT 


45  I F  ( Z-0”’A/> :  6 /  •  So 
50  0l5T2=uI5Tl 
GO  TO  i4u 

60  IF  i  OLZ-U  ( l S'jR  T+ 1 )  )  70 «  ^0 
7w  IF  ( 1  •£-4-b^*vS'j£P)  6c * 140 

80  0l5T2=uI3fl*-2.:>  (OLZ-SOEP)  •'*  ( 5<*\/SDEP)  /  (RSQ  (SK  » \/SDEP)  ) 
GO  TO  14U 

90  0IST2=lI5 ri*2.» (0( I SOPT^l) -50EP)*(V ( ISDPT*1) ♦VSDEP)/ 
1  (R50  (5“  *  -7  ( ISOPT  +  1)  )  ♦R50(5<»VS0EP)  ) 

IF  {  Il2-  T-iSOPT-l )  1  /' <5 ' 1 2 i  L. i, 

100  J=IS0P]+2 

00  110  I=J.ILZPT 
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11:4  i}IST2=”lSr2*£.*C(  I )  /(WSS)(SK»V(I~lI )  ♦RSG(SK*Vf II )  1- 
12*  IF<l*£-4-S**tfLZ?m3«?*l4# 

13*  niST2=='ISi2*2.*C0LZ-GLZ?r»*lS<4VLZ?Tl/*S0CSKtVLZPTl 
CALCULATE  rtAOlUS  AT  TUNING  POINTS 

14*  1-0 1 i ZPT-1 )  1/CVC IZPT-i I -tfZPTI 

KETJ:-<^ 

E^D 


